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Abstract 



The correlation effects in nuclei owing to which the nuclear wave functions are differ- 
ent from the Slater determinants are studied on the basis of the theory of Ref.jlj]. The 
calculated numbers of nucleons out of the nuclear Fermi-surface are in reasonable agree- 



■ ment with the finding from the high-momentum components of the nucleon momentum 

distributions in nuclei. The problems concerning the nuclear binding energy are also 
q ' discussed. 

o 
^ ■ 

1 Introduction 

o 

The nucleon-nucleon interaction still remains one of the central problems of the nuclear structure 
theory. It seems natural to use the free-space forces, but the most of the nuclear structure 
theorists prefer the effective ones instead. The basic motivations are as follows. 

(a) Historically the first one arose from the wide-spread belief in sixties that the free-space 
NN potential has a hard repulsive core. Clearly such interaction does not apply directly because 
of the divergence of finite-order Feynman diagrams. The most popular way to overcome this 
difficulty is to calculate the Brueckner G-matrix and use it as the effective interaction in the 
Hartree-Fock problem. This is the Hartree-Fock-Brueckner approximation, see |2| and the 
references therein for details. But soon it became clear that the description of the two-nucleon 
system, i.e. the deuteron properties and the elastic NN scattering phase shifts below the pion 
production threshold, does not require the hard repulsive core. As a result all the contemporary 
NN potentials are of soft core character, and so the above difficulty does not really exist. In 
such conditions the calculation of the Brueckner G-matrix is not compulsory. 

(b) The renormalization of the free-space interaction due to the medium polarization effects. 
But such effects are treated by conventional methods of quantum many-body theory, and 
therefore the above reason is not the argument for the effective forces. 

(c) The medium QCD-renormalization due to the fact that the quark composition of the 
QCD vacuum is changed in nuclear medium thus leading to the changes of both the mesons 
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and the meson-nucleon vertices 0-0. Such processes can hardly be described in all details 
since the exact theory does not yet exist for the nonperturbative region. Nevertheless, at 
least one exact statement can be done. The QCD-renormalized interaction is the functional of 
nuclear density possessing the following obvious property: it turns into the free-space one in the 
zero-density limit. For this reason it can be represented as the following functional expansion: 

f Q CD-r(riV 2 {p}) = / 2 (|n - r 2 |) + J / 3 (|n - r 2 |, |n - r'|)p(r')dV + 

+ JI / 4 (|n - r 2 |, |n - r'|, In - r"\)p(r')p(r")d 3 r'd 3 r" + • • • (1) 

In this way we conclude that the medium QCD-renormalization is equivalent to the existence 
of many-particle NN forces in addition to the two-particle ones. This conclusion is confirmed 
by the fact that the physics of strong interaction is essentially nonlinear. But the nonlinearity 
automatically leads to many-particle forces. 

The above reasons are grounds for our starting point: both the two-particle and many- 
particle NN forces must be taken into account for the treatment of nuclear structure. 

The three-particle NN forces are indeed included for the calculations of few-nucleon systems 
@. But the following questions arises: is this sufficient for complex nuclei ? This point as 
well as a number of those concerning the nuclear structure might be elucidated if the model- 
independent object would exist in nuclei. 

Our approach is based on the fact that such objects do really exist. As discussed in our 
previous work |]J they are the doorway states for one-nucleon transfer reaction. We showed 
that 

(i) at least the three-particle repulsion and the four-particle attraction must be taken into 
account in addition to the two-particle forces; 

(ii) the nuclear relativity is really existing phenomenon rather than the hypothesis of 
J.D. Walecka @; 

(iii) the dominant contribution to the isovector nuclear potential is provided by the many- 
particle forces; 

(iv) the observed spectra of the doorway states can be used to specify the neutron density 
distributions in nuclei; such specified densities are indeed obtained for the closed-shell nuclei 
40 Ca, 90 Zr and 208 Pb. 

In the present work our approach is applied to the "empirical" studies of nuclear correlation 
effects owing to which the nuclear wave functions are different from the Slater determinants. 
Such effects are treated by a variety of approximate methods since the exact ones do not 
exist. For this reason it is very important to get a model- independent quantitative information 
about the above effects. The possibility provided by our approach is based on the fact that 
the single-particle states of nucleon in nuclear static field (which are just the doorway states 
for one-nucleon transfer reactions JIJ) are correlation-free objects in contrast to the Landau- 
Migdal quasiparticles §| (which include the correlations by definition) and single-particle states 
in nuclear Hartree-Fock calculations with effective forces (where the correlations are implicitly 
included in the phenomenological effective force parameters). So, calculating the correlation- 
free quantities and comparing them with the observed ones we get a quantitative measure of 



2 



the correlation effects. In Sect. 2 this procedure is applied to the nucleon density distributions 
whereas the problems concerning the nuclear binding energy are discussed in Sect. 3. 



2 Density distributions 

The nucleon density distribution in the ground state of nucleus A is (see for the notations) 

(A-l) 

p(r) = (A \^(x)^(x)\A ) = {AoW+WW-lWiA-lWWlAo) = 

j 

{A - 1] r de 

= E = J c ^G(x, X] e). ( 2 ) 

The integration contour C includes the real axis and the infinite radius semicircle in the upper 
part of the complex e plane. As seen from Eq.(l) the density is expressed through the single- 
particle amplitudes of the (A — 1) nuclear states. Expanding them over the complete set of the 
doorway states, 

= Ec] A Va(x) (3) 

A 

and putting Eq.(3) into Eq.(2) we get 

P(r) = E ™a|^a(x)| 2 + 2]T pxui>i(x)i) v (x) , (4) 



where 

(A-l) (A-l) (A-l) 

n, = Pxx = E l^ A) ! 2 = E # } > PX»= E Ci W 'cf (5) 



J J 3 



(actually the coefficients Cj are real quantities if the parity violation due to weak interaction 
is disregarded). The diagonal elements p\\ = n x are the doorway state occupation numbers. 
Indeed, the particle number is 



N 



J p(r)d 3 r = E n A , (6) 



since the nondiagonal elements do not contribute because of the orthogonality between ip\ and 
ip u . The quantities n\ and p\ v obey the following limitations: 

< n x < 1 , \p Xu \ < ^(n x + n u ). (7) 

The first follows from the fact that the doorway states are distributed over the actual ones of 
both the (A — l) and (A + l) nuclei, see Eq.(13c) in Ref. ffl, whereas the second is a consequence 
of the Cauchy-Bunjakowsky inequality. 

The facts n x < 1 and p Xv ^ reflect the Fermi surface smearing due to the correlation 
effects. The quantities n x and p Xv carrying the quantitative information about these effects 
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Figure 1: The observed (full line) and correlation-free (dashed one) proton density distributions 
in 40 Ca. 

can be found "empirically" by considering the relation (4) together with the limitations (7) as 
equations for n A and p\ v . The results will be published elsewhere. 

As a result of the correlations a part of nucleons is out of the nuclear Fermi-surface. The 
number N out of such nucleons is calculated by comparing the observed density distribution with 
the correlation-free one 

Pcf(r) = £0(^-£ A )hM*)r. (8) 

A 

Bp is the Fermi level energy. Such comparison is shown in Fig.l for the proton density distri- 
bution in 40 Ca. As seen from the figure the correlation-free density contains more nucleons in 
the inner region < r < than the observed one, Tj is the intersection point. The situation 
in the outer region r,i < r < oo is clearly opposite since both densities correspond to the same 
number of nucleons. The number of redistributed nucleons is 

n oo 
N ut = 4tt J \p c f{r) - p(r)]r 2 dr = 4vr J [p(r) - p cf (p)]r 2 dr . (9) 
o n 

Ti is the intersection point: p c /(?"i) = piji). This is just the N out because the only reason for 
the redistribution is the Fermi-surface smearing due to the correlations. 

The N out numbers in doubly closed-shell nuclei are shown in Table 1. As seen from the 
table these nucleons are a rather appreciable part of the total mass number. To our knowledge 
this fact was first mentioned by Frankfurt and Strikman [JTOj on the basis of the analysis 
of high-momentum components (i.e. those with k > 300MeV/c) of the nucleon momentum 



distributions in nuclei. According to their most recent results for these data |TI[ the A out /A 
ratio is (20 ±3)% for heavy nuclei. This is in reasonable agreement with our results. 

It is worth mentioning that our calculations do not tell anything about the nature of the 
underlying correlations whereas the high-momentum tales of the momentum distributions arise 
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Table 1: iV oll4 numbers in doubly closed-shell nuclei 
16 4U Ca sn & m W 

N out 1.10 3.15 7.34 13.15 

Z out 1.15 3.53 6.46 13.22 

A out 2.25 6.68 13.80 26.37 

A out /A, % 14 16.7 15.3 12.7 



from the NN interactions at short distances [JO], |TT[. Therefore the reasonable agreement 
between the two results gives rise to the conclusion that the main reason for the Fermi-surface 
smearing in doubly closed-shell nuclei is due to the short-range correlations. 



3 Binding energy 



We also can calculate the part of the nuclear binding energy which is caused by the motion of 
nucleons in nuclear static field. Such static energy partly includes the correlation effects since it 
is expressed through the observed nucleon density distributions. The comparison between the 
observed binding energy Sb and the static one E st gives the measure of the proper correlation 
energy of nucleus. 

To make this point clear let us derive the exact expression for the binding energy. Following 
the procedure of Ref . |§ we get 

S b = JdxJ -^trk x G(x,x;e) + ^JJdxdxx J J 7^~p -M|r - r x |); uj)K 2 {x, xi; x, xx] e, £x) + 



+ 7T 



+ j 



x K. 



7^ JJJ dxdxidx 2 JJJ £ (2ri)3 2 f 3 (\ r ~ ri l' l r ~ r2 \i u > Ui)K 3 (x,x 1 ,x 2 ; x, xi, x 2 ; e, s 1 ,s 2 ) + 

c 

\ffff dxdx x dx 2 dx z JJJ J ^^^|^/ 4 (|r - r^, |r - r 2 |, |r - r 3 |; u, io u u 2 ) x 

c 

\(x,xi,x 2 ,x 3 ;x,x 1 ,x 2 ,x 3 ;e,e 1 ,£ 2 ,£ 3 ) , (10) 



where K(xx, . . . x n ; x[ . . . x' n ; E\ . . . e n ) are the n-particle Green functions. We accounted for 
the three- and four-particle forces in addition to the two-particle ones as well as the possible 
dependence of the interactions on the appropriate energy transfers Wj. In these terms the 
single-particle Green function is 

(e - k x )G(x,x';e) = 5(x - x) + J dx% J 7^/2 (l r ~ r ih u)K 2 {x, x x \ x', x[; e, + 



+ 



JJdx x dx 2 J J /3(|r-ri|,|r-r 2 |;a;,a;i)lif3(a:,xi,X2;a/,zi,X2;e,e2) + 
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+ JJJ dx x dx 2 dx 3 JJJ £l , 2 £ L £3 /4(lr - ril, [r-r 2 |, |r- r 3 |;a>,a>i,a> 2 ) x 

c 

x K 4 (x,x 1 ,x 2 ,x 3 ;x',xi,x 2 ,x 3 ;e,e 1 ,e 2 ,e 3 ) . (11) 
On the other hand the Dyson equation is 

(e — k x )G(x, x'\ e) = 5(x — a/) + J dx 1 M(x,x 1 ;e)G(xi,x';e) , (12) 

and therefore 

J M(x,x 1 ;e)G(x 1 ,x';e)dx 1 = J dx x J ^f 2 (\r -r 1 );uj)K 2 (x,x 1 ;x\xi;e,ei) + 

c 

+ /7 dx x dx 2 JJ d (^ d \l M\ T ~ r iU r ~ r 2[; ^ k>i)if 3 (:E, x 1; z 2 ; x', x x , x 2 ; e, e h e 2 ) + 

c 

+ JJJ dx x dx 2 dx 3 JJJ d£ ^ £2< ^ 3 U(\r - ri|, |r - r 2 [, [r - r 3 |; a>, a>i, u 2 ) x 

c 

X i^ 4 (a;,a;i,X2,X3;x / ,xi,X2,X3;£:,ei,e2,£3) • (13) 
As seen from Eqs. (10) and (13) the binding energy can be written as 

£ b = J J dxdx' J (trk x 5{x — x') + M(x, x'\ sfj G(x r , x; e) — 

c 

\l I dxdXl IS (2 £ )2 ^ r ~ T ^ U ) K ^ x^x^^e, £i) - 

c 

2 [ffj j j /"/"/" dede\de 2 



ujdxdxxdx^jj / 2 ^ 3 /aflr - r x |, |r - r 2 |; a;, ^i)if 3 (x, g 1; x 2 ; x, g l8 x 2 ; £i, e 2 ) - 

c 

^ JJJJ dxdxtdx 2 dx 3 JJJJ ^^^^/ 4 (|r - r x | , |r - r 2 |, |r - r 3 \,u,Ui,(v 2 ) x 

c 

x K 4 (a;,a:i,x 2 ,a; 3 ;a;,xi,a; 2 ,a; 3 ;£,£i,£ 2 ,£ 3 ) . (14) 

Taking into account the spectral representation of G(x,x';e), see Eq.(5) in Ref.0], and per- 
forming the integration over e we get the following expression for the first term in the rhs: 

(A-i) , . 

e i . = E dxty+(x) (k^jix) + Mixix'^E^jix^dx') . (15) 

i 

As follows from the Dyson equation (12) and the above spectral representation the amplitudes 
*&j(x) obey the equation 

k x tyj{x) + J M{x,x';Ej)y!j{x')dx' = E^^x) (16) 

and therefore 

4 1] = E e j J \*j(x)\ 2 dx = EE M A) ( 17 ) 

3 A j 



(we accounted for the expansion (3)). 

In general case the many-particle Green functions K n (xi, . . . x n ; x[, . . . x' n ; . . . e n ) obey 
the infinite system of integro-differential equations the Eq.(ll) being the first one. As men- 
tioned above the exact methods of solution do not exist. It should be also mentioned that the 
approximate methods, see [EJ and the references therein, are developed for the instantaneous 
two-particle forces only. For these reasons the exact calculation of nuclear binding energy is 
impossible at present. Therefore the calculation of the static energy is of importance. 

Both the static nuclear field and the static energy are obtained from the above relations by 
putting 

n 

-Kn{%li ■ ■ ■ %n'i %li ■ ■ ■ %n'-> ^1) • • • ^n) = W G(X{, Xf, Si) , (18) 

i=l 

i.e. neglecting the difference between the many-particle Green functions and the unsymmetrized 
products of the single-particle ones. This difference arises from both the antisymmetrization 
and the higher-order terms of the perturbation theory, i.e. just the effects leading to the proper 
correlation energy. Putting Eq.(18) into Eqs.(13) and (14) we get (all the energy transfers Ui 
vanish in this case) 

U s t{r) = f / 2 (|r-ri|)p(ri)d 3 ri+ If / 3 (|r - n|, |r - r 2 |)p(ri)p(r 2 )d 3 rid 3 r 2 + 

+ //// 4 (|r-r 1 |,|r-r 2 |,|r-r 3 |) P ( ri ) P (r 2 Wr3)A 1 A 2 A3 ; (19) 

{A - 1] r ^ 1 rr 

£st= E / ^ + (x)(^ + ^(r))^(x)rfx--yy/ 2 (|r-r 1 |)p(r)p(r 1 ) C / 3 r C / 3 r 1 - 

3 

\fff ' f 3 (\r-r 1 \,\r-r 2 \)p(r)p(r 1 )p(r 2 )d 3 rd 3 r 1 d i r 2 - (20) 

I IIII ~ ri '' ' r ~ r2 '' ' r ~ r 3\)p( r )p( r i)p(r 2 )p(r 3 )d 3 rd i r 1 d 3 r 2 dh 3 . 
As follows from the relations (3) and (5) the first term of the rhs is 

& = J2n x e x . (21) 

A 

But as seen from Eq.(17) and the sum rule 

(A-l) (A+l) 

^ = E E A X) + E M x) (22) 

j k 

for the doorway state energies, see Eq.(14c) of Ref.[Q, the quantity contains unphysical 
contributions from the energies of the (A+l) nuclear states. To eliminate them let us divide 
into the two parts 

e$ = e< + s> , e<= E ^a^a , e> = E ^a^ a , (23) 

A<F A>F 

including summations over the states with S\ < Ef and e\ > Ep respectively. The doorway 
states entering £^ are mainly distributed over the states of the (A — 1) nucleus the first term 
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of Eq.(22) thus giving the dominant contribution to the E\ values in this case. The small 
unphysical contribution from the second term of (22) arises from the states of the (A + 1) 
nucleus with the same quantum numbers as those of the low- lying states of the (A — 1) nucleus. 
Such states lie either in the continuum or near its border. Therefore the reasonable estimate 
for their energies is -E^ ~ the unphysical contribution thus being negligible for this case. 

The situation is opposite for £f t because the doorway states with e\ > Ef are mainly 
distributed over the states of the (-4+1) nucleus, the dominant contribution to the E\ values thus 
being unphysical from the viewpoint of the binding energy. In such conditions it is reasonable 
to use a more appropriate relation for £, 



(A-l) 



£ > 

° St 



E E 



(24) 



\>F 



So the (A — 1) nuclear states with the same quantum numbers as those of the low-lying states 
of the (A + 1) nucleus contribute in this case. There are no such states among the low- lying 
ones of the (A — 1) nucleus, and therefore 



E S <E F = £ (A)-£ g (A-l) , 
see Ref. |I] for the details. Nevertheless it is reasonable to use the estimate 

Ej = E F , 

providing the least absolute value for the £ ^ quantity. In this way we get 



(25) 



(26) 



(A-l) 



£> 

° St 



e f j: e 



\>F j 



Ef^tix = EpNout 

\>F 



(27) 



Accounting for the fact that the many-particle forces are introduced as the contact ones in the 
model-independent approach of Ref. M , we finally get (see Ref. M for the notations) 



£ «t = E E ™a£a + N out E F 

i=n \\<F 

1 f [S(r)p s (r) + V u {r)p{r) + S~{r)pj(r) + V p {r)p-{r) + C{r)p ch {r)) dh 



2 3 1 r -i 

-a 3 p 3 (r) + -a 4 p 4 (r) + - [% p{r) + a^p 2 (r)\ (p"(r)) 2 



We used the following ansatz for the occupation numbers 



1 - 



( £A -/i) 2 + C 2 



E^A 
A 



N. 



E ™A 

A>F 



N. 



out > 



(28) 



(29) 



putting N out = N out , see Sect. 2. This is incorrect because the sum J2\>F n \ is different from 
N out , Eq.(9). Indeed, both diagonal and nondiagonal elements of the density matrix, Eq.(4), 
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Table 2: The static and observed binding energies in MeV together with the dominant contri- 
butions to £ st . Those from the Coulomb and the isovector terms are not shown, but they are 
included into the E st value. The notations correspond to the terms of Eq.(28). 





°st 


2 ^ 


-\v u 


— |d3 




Est 


s b 


16Q 


-345.6 


+ 1528.4 


-1224.0 


-451.2 


+461.8 


-47.0 


-127.6 


40 Ca 


-991.7 


+4456.0 


-3570.8 


-1346.0 


+1411.3 


-120.7 


-342.0 


90 Zr 


-2251.6 


+10699.0 


-8565.6 


-3236.5 


+3386.0 


-234.2 


-783.9 


208p b 


-4998.2 


+26463.4 


-21255.6 


-8190.9 


+8694.0 


-233.3 


-1636.5 



contribute to the N out values. This shortcoming will be corrected in future by calculating the n\ 
and pxu values, see the discussion in Sect. 2 (of course, the ansatz (29) will become unnecessary). 

The results of the calculations are shown in Table 2. As seen from the table the static 
energy is the sum of a number of contributions with different sign, the dominant ones greatly 
exceeding the total static energy as well as the observed binding one. This is the source of the 
ambiguities because all disregarded effects may be of importance in such a situation. 

One of such effects is the finite range of the many-particle forces. The possible way of 
the estimation is illustrated for the three-particle forces. As seen from Eq.(10) the original 
contribution is 

£3 = ~ III /3(|r - ri|, |r - r 2 \)p(r)p(r 1 )p(r 2 )d 3 rd 3 r 1 d 3 r 2 
p{r)m, vM\* + t\)p(\r + r]\)d 3 rd 3 $d 3 r, - 



p(r) + £Vp(r) + -iiiAd kP {r) + ■ 



p(r) + 77 Vp(r) + -rijriedjdepir) + 



d 3 rd 3 £d 3 ri 



(30) 



p{r) + -£ 2 Ap(r) + 
6 



p{r) + -r] 2 Ap(r) + ■ 
6 



-a 3 I p(r)[p(r) + r^Ap(r) + • • fd s r 



d 3 rd 3 £d 3 r) 

If 2 r 

-a 3 J p 3 {r)d 3 r+ -a 3 J p 2 (r)Ap(r)d 3 r 



where (see Eq.(37) of Ref.0) 



«3 



^Hi^dHd 3 ^ . (31) 



In this way we get 



5 r? S st = r 2 J (Ja 3 p 2 M + ^a 4 p 3 (r)^) Ap(r)d 3 r = -Dr, 
D = J (la 3 p(r) + ^a 4 p 2 (r))(p'(r)) 2 rf 3 r. 



(32) 



The calculated D values are shown in Table 3. The reasonable value of the range parameter 
is the ui meson Compton wavelength, i.e. a typical scale of the strong interaction. So r 2 = 
m~ 2 = 0.1 fm 2 , the effect thus being negligibly small. 
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Table 3: The many-particle finite range effect 



167 



O 40 Ca 9U Zr 2U8 Pb 



D, MeV-fm" 2 -12.12 -3.29 22.02 34.0J 



Table 4: The contributions to the static field 



u 2 


u 3 




Bonn B -83 


+96.5 


-104 


OSBEP -82 


+97 


-107 



The second source of the possible ambiguity is seen from the results of Ref.0 for the 
contributions of two-particle, three-particle and four-particle forces to the nuclear static field 
which are shown in Table 4, see Eqs. (44) and (45) of Ref. |l|. 

As follows from the table the convergence of this sequence is not seen. This observation sug- 
gests the possible existence of the contributions from higher many-particle forces, and therefore 
the whole sequence should be summed up. But this can be done only within some reasonable 
model for the many-particle forces including all higher-order ones as well as the finite range 
and the mechanism for the saturation of nuclear density. The corresponding investigation is in 
progress. 
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